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Section  1 


INTRODUCTION 


Given  signals  from  two  spatially  separated  hydro 
phones  which  receive  energy  from  the  same  broadband  source 
it  is  possible  to  measure  the  travel-time  difference  from 
the  source  to  those  phones.  The  success  in  doing  such 
a measurement  is  dependent  on  several  practical  factors: 

• the  extent  to  which  the  signal  on  the  two 
phones  is  spatially  coherent,  and 

• the  bandwidth  of  the  signal. 

For  the  problem  of  interest  here,  the  source  will  be  with- 
in direct  path  range  of  the  two  receivers  (which  are  no 
more  than  a few  thousand  feet  apart)  and  so  it  will  be 
assumed  that  the  same  signal  is  received  on  both  phones 
and  lack  of  coherence  is  not  a problem.  The  significance 
of  bandwidth  will  be  discussed  below. 

The  problem  of  interest,  alluded  to  above,  is 
that  of  establishing  the  travel-time  difference  as  a 
function  of  time  for  a surface  ship  transiting  past  bot- 
tomed hydrophones.  A formalism  is  developed  in  the  next 
section  which  permits  establishing  this  time  difference 
from  an  analysis  of  the  normalized  cross  covariance  of 
the  data  from  pairs  of  phones. 


1-1 


Section  2 


RELATIONSHIP  OF  TRAVEL-TIME 
DIFFERENCES  TO  CROSS  COVARIANCE 

2.1  Some  Basic  Relationships 

We  begin  by  expressing  the  signal  as  received 
on  a single  phone  from  some  number,  M,  of  discrete  sources 
as 

M 

x(t)  = 5]  a. (t)  (1) 

i=l  ^ 

^ ll 

where  a^(t)  is  the  signal  from  the  i—  source.  The  sig- 
nal received  by  a second  phone  some  distance  from  the 
first  is,  then, 

M 

y(t)  » X)  a (t-T.)  (2) 

i-1  ^ 

where  the  signal  from  the  i—  source  is  delayed  by  some 
time  from  its  arrival  at  the  first  phone.  Several 
assumptions  are  implicit  in  Eqs . (1)  and  (2): 


• there  is  perfect  coherence  between  the  two 
phones ; 


• all  the  sources  of  significant  energy  are 
discrete  (we  are  interested  in  that  part 
of  the  spectrum  which  is  dominated  by 
broadband  noise  from  surface  ships). 

It  will  further  be  assumed  that  signals  from  two  differ- 
ent sources  are  uncorrelated,  which  can  be  expressed 
here  as 


<a^(t)aj(t+A)>  = 6^jR^(A) 


(3) 


where  R.(A)  is  the  autocorrelation  of  the  signal  from 

1“  h ^ 

the  i—  source  and  the  brackets  denote  statistical  expec- 
tation . 


It  follows  from  Eqs . (1),  (2)  and  (3)  that 
the  autocorrelation  for  x(t)  and  y(t)  and  the  cross 
correlation  between  them  are  (assximing  the  processes  are 
stationary) 

R^(A)  = <x(t)x(t+A)> 


= 53  <a.  (t)a  .(t-t-A)> 

i,  J ^ 


-LRi(A)  (4) 
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Ry(A)  = <y(t)y(t+A)> 


= H a^(t-T^)aj(t-T^+A) 


= ERi(M 

i 


= R (A) 

X 


(5) 


Rjjy(A)  = <x(  t)y(t+A)> 


= E 

ij 


<a.  ( t)a'.(  t-T  .+A)> 
J 3 


= £ r.(a-t.). 
. 1 
1 


(6) 


Because  the  signals,  x and  y,  are  acoustic  pressures, 
they  have  zero  mean  and  so  the  normalized  cross  covariance 
can  be  written  directly  as^ 


R (A) 

p (A)  = j 

(R^(o)R  (o))* 


E 


E 

i 


Ri( A-ti) 

R^Co) 


(7) 
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It  will  further  be  assumed  that  the  processes,  x and  y, 
are  ergodic  so  that  the  various  correlation  functions 
can  be  estimated  from  a sample  time  series: 

■T 

Rjj(A)  - ^ / xCt)x(t+A)dt  (8) 

o 

T 

R (A)  - h f y(t)y(t+A)dt  (9) 

y o 


(A) 


x( t )y( t+A)dt . 


(10) 


For  the  case  where  A<<T,  it  follows  from  Eqs.  (8),  (9) 
and  (10)  that  the  correlation  functions  also  can  be 
estimated^  from  the  Fourier  transform  of  either  a power 
density  spectrum  or  a cross  spectral  density: 

R^(A)  = F"^  [lS^(a))|2]  (11) 

(12) 

(13) 

where 

Sj^(ui)  = F[x(t)] 

i 

i 
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Ry(A)  = F"^[|Sy(a))|2] 

R^y(A)  = F"^[S^(a))Sy(u))*] 


I 


* Sy{co)  i F[y(t)] 


F,  F ^ and  * denote  the  Fourier  transform,  the  inverse 
t Fourier  transform  and  complex  conjugation,  respectively. 

Finally,  note  that  the  auto  correlation  functions  R (A), 
Ry(A)  and  I^^(A)  when  evaluated  for  A = 0 are  estimates 
of  the  total  intensity  for  the  signals  x(t),  y(t)  and 
a^(t),  respectively. 


i 

I 


I 


i 

1 


I 


For  a broadband  signal  with  an  essentially 
flat  spectriim  of  width  the  auto  correlation  function 
R.(A)  will  be  a function  which  peaks  at  A = 0 and  has 
a nominal  width  of  2ir/Q.  So  the  expression  ^ R^(A-t^), 
from  the  numerator  of  Eq.  (7),  might  appear  as  shown  in 
Fig.  1.  In  this  stylized  picture,  it  has  been  assumed 
that  the  widths  of  the  R^  functions  are  comparable  and 
the  time  delays  are  such  that  most  of  these  peaks  are 
separable.  If  the  i—  peak  is,  indeed,  distinguishable 
from  the  rest,  then  Eq.  (7)  evaluated  at  A = is,  approxi- 
mately. 


R.(0) 

ER.(O) 


(14) 


from  which  it  is  clear  that  at  the  time  delay  x^,  the 

cross  covariance  is  just  the  ratio  of  the  intensity  of 
t tl 

the  i—  signal  to  the  total  intensity  on  either  of  the 
phones.  It  is  intuitively  clear  at  this  point  that  for 
the  case  of  one  nearby  broadband  source  and  many  at  sub- 
stantially greater  distances  (and  varied  time  delays). 
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the  nearby  source  should  show  up  as  a distinct  peak  in 
the  covariance  at  a time  delay  which  can  be  measured 
with  some  accuracy  equal  to  or  less  than  In 

section  2.3  below  an  estimate  of  what  could  be  expected 
is  made  for  a typical  case  of  interest. 

2.2  Compensation  for  Doppler  Shift  and  Discrete  Lines 

Because"  of  the  short  ranges  involved  for  ships 
within  direct  path  range  of  the  phones,  it  is  likely 
that  a relative  Doppler  shift  will  exist  which  must  be 
compensated  before  computing  the  cross  correlation  of 
the  two  phones.  Because  the  cross  correlation  is  com- 
puted from  the  cross  spectral  density,  it  is  convenient 
tp  adjust  for  the  relative  Doppler  shift  in  the  frequency 
domain.  Fig.  2 illustrates  the  way  this  shift  has  been 
implemented.  Since  the  data  are  processed  in  digital 
form,  one  has  N sample  points  of  each  spectrum  which 
is  bandpass  filtered  as  shown.  This  spectrum  is  resam- 
pled by  linear  interpolation  of  the  amplitude  and  phase 
to  obtain  N samples  on  the  frequency  shifted  axis  where 

CO  = ( 1+6  )co  ^ . 

The  unshifted  frequency  is  co,  the  shifted  frequency  is 
co^  and  6 is  the  fractional  shift  in  frequency.  To  obtain 
the  necessary  out-of-band  values  along  the  co'  axis  which 
extend  beyond  the  last  value  on  the  co  axis,  the  unshifted 
spectrum  is  zero-filled  as  necessary.  By  shifting  either 
S^(co)  or  Sy((o),  the  two  possible  relative  Doppler  shifts 
(up  or  down)  can  be  compensated. 
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The  actual  spectrum  of  a ship  wil_  look  as 
shown  in  Fig.  3,  containing  both  a broadband  background 
and  discrete  lines.  If  a significant  fraction  of  the 
total  intensity  is  in  the  lines,  then  they  can  compli- 
cate the  cross  correlation  function  by  the  addition  of 
discrete  sinusoidal  components,  one  for  each  line  in 
the  spectrum  and  having  the  frequency  of  that  line. 
Accordingly,  the  formalism  for  constructing  the  cross 
correlation,  R^y(A),  includes  the  capability  to  elimi- 
nate any  strong  line  components  from  S (co)  and  S (w) 

X y 

prior  to  computing  the  cross  correlation.  This  is  accom- 
plished by  a sliding  window  which  compares  the  amplitude 
value  at  its  midpoint  with  the  median  value  for  the  win- 
dow (exclusive  of  the  point  being  compared);  if  the 
amplitude  of  the  tested  point  exceeds  the  median  ampli- 
tude by  more  than  a specified  amount,  xt  is  replaced 
by  the  median  amplitude. 

2.3  An  Estimate  of  What  to  Expect  for  a Case 

of  Interest 

Consider  a case  where  the  ambient  spectral 
background  has  a level  of  80  dByPa/Hz  in  the  absence  of 
nearby  ships.  Given  a nearby  ship  which  has  a radiated 
spectral  level  of  160  dByPa/Hz,  how  far  away  (assuming 
spherical  spreading  along  the  direct  path)  could  this 
ship  be  and  still  yield  a peak  in  the  normalized  cross 
spectral  density  of  at  least  0.4?  Assuming  Eq.  14  to 
be  an  appropriate  approximation  for  this  case,  we  re- 
quire 


(15) 
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where 


10  log  I = 160  - TL 
10  log  N = 80 

and  we  want  TL  such  that  Eq . (15)  is  satisfied.  This 
implies 

TL  ^81.4 

or,  assuming  spherical  spreading,  that  the  range  satis- 
fies 

r ^ 11  km 

Since  the  direct  path  range  for  a bottomed  phone  at  a 
depth  of  2 - 4 km  is  generally  less  than  11  km,  this 
analysis  suggests  that  any  ship  within  direct  path  will 
result  in  an  obvious  peak  in  the  cross  covariance  under 
the  representative  conditions  given  above. 
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Section  3 


SUMMARY 


As  a step  in  tracking  surface  ships,  a forma- 
lism has  been  developed  to  measure  the  time  history  of 
travel-time  differences  between  pairs  of  phones  by  cross 
correlating  on  the  broadband  output  of  ships  within 
direct  path  range  of  the  phones.  Provisions  have  been 
made  for  compensating  relative  Doppler  shift  and  elimi- 
nating (when  necessary)  strong  line  components  orior 
to  computing  the  cross  correlation.  By  assuming  typical 
values  for  the  ambient  noise  and  the  ship  radiated  level, 
it  was  shown  that  a ship  within  direct  path  range  could 
be  expected  to  yield  a peak  in  the  normalized  cross 
covariance  of  at  least  0.4.  These  algorithms  have  been 
incorporated  in  a computer  code,  as  described  in  Reference  2. 

\ 
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ance  function 


Figure  2.  Relationship  of  shifted  to  unshifted  spectrum  levels 
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